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We provide a new derivation of the anisotropics of the cosmic microwave background (CMB), 
and find an exact expression that can be readily expanded perturbatively. Close attention is paid 
to gauge issues, with the motivation to examine the effect of super-Hubble modes on the CMB. 
We calculate a transfer function that encodes the behaviour of the dipole, and examine its long- 
wavelength behaviour. We show that contributions to the dipole from adiabatic super-Hubble modes 
are strongly suppressed, even in the presence of a cosmological constant, contrary to claims in the 
literature. We also introduce a naturally defined CMB monopole, which exhibits closely analogous 
long-wavelength behaviour. We discuss the geometrical origin of this super-Hubble suppression, 
pointing out that it is a simple reflection of adiabaticity, and hence argue that it will occur regardless 
of the matter content. 

PACS numbers: 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 

The anisotropics in the cosmic microwave background 
(CMB) reveal a great deal about our Universe, since they 
persist essentially unscathed from the epoch when fluc- 
tuations were well described by simple linear theory. The 
comparison of CMB observations with theory has become 
a mature subject, and has played an important role in 
forming our current understanding of the Universe (see, 
e.g., [1]). 

Critical to that comparison is the accurate theoretical 
calculation of the anisotropies. Since the pioneering work 
of Sachs and Wolfe [2] , the theoretical anisotropies have 
been refined and recalculated using different formalisms 
many times (see, e.g., [3, 4, 5, 6, 7, 8, 9]). Accurate calcu- 
lations are now readily available via public code packages 
such as CAMB [10, 11]. 

In the present work, we revisit the calculation of 
anisotropics. While our results may not lead to more 
accurate or efficient calculations, wc hope that they will 
help to clarify some of the conceptual issues surround- 
ing the calculations. In particular, our approach makes 
explicit the physical meaning of the various contribu- 
tions to the anisotropies. Crucial to this is our use of 
the covariant approach to cosmology (see [12, 13] for re- 
views), which is ideal for writing exact solutions and for 
physical clarity. We present a remarkably simple but ex- 
act expression for the anisotropy, which applies to arbi- 
trary spacctimcs and includes the effects of tensor as well 
as scalar perturbations and any line-of-sight integrated 
Sachs- Wolfe (ISW) effect. This general result can be 
readily expanded perturbatively, and here we turn to the 
metric formalism for computational efficiency and show 
that we recover previous results in the literature. 

A main motivation for our work is in examining 
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the behaviour of the anisotropies due to super-Hubble 
fluctuations (we use the terms "super-Hubble", "long- 
wavelength" , and "large-scale" interchangeably, to mean 
scales larger than the current Hubble or last scatter- 
ing radius), where gauge issues are paramount. This 
question has been examined before in the context of 
the Grishchuck-Zel'dovich effect [14], which describes the 
large-angular-scale anisotropics that result from super- 
Hubble modes. In the context of a matter-dominated 
universe with adiabatic fluctuations, it was shown that 
the CMB dipole receives strongly suppressed contribu- 
tions from long-wavelength modes. A claim was made in 
Ref. [15] that this suppression would not occur in mod- 
els with cosmological constant, so that we could "see" 
very long-wavelength structure in the dipole. It was also 
found that in the presence of isocurvature perturbations, 
the suppression may not occur (see [16] and references 
therein) . 

To study this issue, we construct a transfer function 
that describes the scale dependence of contributions to 
the dipole. Working by analogy, wc carefully define a 
CMB monopole perturbation, and find its transfer func- 
tion. As is well known, such a monopole cannot be ob- 
servable, but we show that its variance is well defined the- 
oretically. Our definitions have simple interpretations: 
the dipole measures the departure of radiation and mat- 
ter comoving worldlines, while the monopole measures 
how well radiation and matter constant-density hyper- 
surfaces coincide. The usefulness of the monopole will 
be in examining its long-wavelength behaviour, where it 
will help to clarify the dipole case. 

We show that the contributions to both dipole and 
monopole vanish for large scale sources, even in the pres- 
ence of a cosmological constant. We close by pointing 
out that this is a direct consequence of adiabaticity, and 
hence that this result is expected to hold regardless of 
the matter content of the Universe, unless isocurvature 
modes are present. 

Another potential reason that a careful treatment of 
the monopole may be of interest involves the measure- 
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ment of the mean CMB temperature and its relation to 
constraints on other cosmological parameters. The mean 
temperature is currently measured to a precision of a few 
parts in 10 4 [17]. It has been pointed out that the pre- 
cision of this measurement could be improved by nearly 
two orders of magnitude with currently available tech- 
nology [18]. Such a measurement would reach the naive 
cosmic variance limit of a part in 10 5 as suggested by the 
observed amplitude of fluctuations (see [19] for a related 
discussion). It would then become necessary to be very 
careful about exactly what information the mean tem- 
perature measurement is giving us, and about the na- 
ture of monopole fluctuations. Relevantly, recent studies 
have examined the importance of the mean temperature 
measurement to our ability to constrain the cosmological 
parameters [20, 21]. 

Since the covariant approach to cosmology is essen- 
tial to this work, we begin in Sec. II with a summary of 
the required formalism. Next, in Sec. Ill we present the 
derivation of the Sachs- Wolfe effect, beginning with an 
exact result before specializing to first order and recov- 
ering previous results. In Sec. IV we present calculations 
of the dipole and monopole transfer functions, and we 
examine their long-wavelength behaviour in Sec. V. Fi- 
nally we discuss our results in Sec. VI. The Appendices 
summarize relevant material in the metric formalism, and 
demonstrate both the gauge invariancc and the gauge de- 
pendence of our results. We use signature (—,+,+, +), 
and greek indices indicate four-tensors, while latin indices 
indicate spatial three-tensors. 



II. COVARIANT COSMOLOGY 

This section will contain a brief summary of the ele- 
ments of the covariant approach to cosmology (see, e.g., 
the reviews [12, 13]) that will be needed in the following 
sections. (A collection of results from the metric-based 
approach to cosmological perturbations, which will also 
be needed, is presented in Appendix A.) Fundamental to 
the covariant approach is the notion of a congruence of 
worldlines, also known as a threading of the spacetime or 
sometimes as a choice of frame. This is a family of world- 
lines such that exactly one worldlinc passes through each 
event. A useful example is the congruence of comoving 
worldlines, when it is well defined. A timclikc congruence 
is described by a vector field u p , tangent everywhere to 
the worldlines. We will assume that is normalized, 
u^u^ = — 1. At each event we can define the spatial 



projection tensor 
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which projects orthogonal to u^. Hypersurfaces orthog- 
onal everywhere to u p will exist if the twist of the con- 
gruence (defined below) vanishes [22], in which case h^ v 
is the (Riemannian) metric tensor for those spatial hy- 
persurfaces. 



It is useful to describe the geometrical properties of the 
congruence by the covariant derivative u KV . By virtue of 
the normalization condition, this derivative satisfies 



0. 



(2) 



The derivative can be decomposed into parts parallel to 
and orthogonal to in its second index using Eq. (1), 
giving 



(3) 



The temporal part can be written in terms of the accel- 
eration of the worldlines, defined by 



(4) 



The field a M measures the departure of the worldlines 
from geodesic. The spatial part h p v ii„. p can be decom- 
posed into trace, symmetric trace- free, and antisymmet- 
ric parts via 



= h p ^u„. r 



er„„ = h p , u x 0h„ v , 
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(5) 

(6) 
(7) 



The scalar 9 and tensors and w^j, measure the local 
rates of expansion, shear, and twist of the congruence, 
respectively. Combining Eqs. (3) to (7) we have 

u^v = -dh^ + (J + uj^ - a^Uy. (8) 

For the case of a homogeneous and isotropic Friedmann- 
Robertson- Walker (FRW) cosmology, if we choose to 
be comoving then = lo^ = a M = and 9 = 3H, 
where H = a/a is the Hubble rate, and a the scale factor. 

We will also need two derivatives. For an arbitrary 
tensor X , the covariant time derivative is defined by X = 
A ;AI it M , and gives the proper time derivative along u^. 
The symbol represents the spatial (orthogonal to u p ) 
covariant derivative defined using the spatial metric h^. 
For example, 



D..X" 



h x h 1 * X p 



(9) 



for any tensor X^ orthogonal to u^. 

A congruence u p can be used to describe the matter 
content as observed locally by a family of observers, given 
the energy-momentum tensor T^. The energy density 
p, momentum density g„, pressure P, and anisotropic 
(shear) stress tt^ , as viewed locally by an observer with 
four- velocity u p , are defined by the projections 



T^vfu" = p, 
T^h\ = -q K , 

T lu/ h\h v lt =Ph XH +ir Xlt , 
where -Kn V is defined to be trace-free. 



(10) 

(11) 
(12) 
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III. SACHS- WOLFE EFFECT 

In this section, we will provide a calculation of the 
anisotropics of the CMB radiation. The main goal will 
be to clarify the issues of gauge and frame choice, which 
will be critical to properly describing the dipole and 
monopole anisotropics in the next section. For this rea- 
son, it will not be necessary to consider here the finite 
thickness of the last scattering surface, or the effect of 
reionization on the anisotropics, which have negligible 
effect on the largest-scale anisotropics. Thus we assume 
tight coupling in the baryon-photon plasma before last 
scattering, followed by instantaneous recombination of 
matter and free streaming (i.e. unattenuated geodesic 
evolution) of the radiation. Apart from this approxima- 
tion the calculation will be exact. 



A. Exact expression 

We are making the approximation that the CMB ra- 
diation is emitted abruptly when the local plasma tem- 
perature drops below some value Te at which recombi- 
nation occurs (E for "emission"), and thereafter trav- 
els freely. This temperature is defined with respect to 
the frame (or congruence) which is comoving with 
the plasma (i.e. for which the plasma momentum den- 
sity vanishes). The spacclike hypcrsurfacc defined by the 
moment of recombination will be called the last scatter- 
ing hypcrsurfacc, Els, while the intersection of Els with 
an observer's past light cone defines the two-sphere com- 
monly called the observer's last scattering surface (LSS). 
Via the Stefan-Boltzmann law, Els must be a hypersur- 
face of constant radiation energy density Pm, with the 
density again defined with respect to the comoving con- 
gruence u M [23] . For adiabatic perturbations at last scat- 
tering and on large scales, Els will also be at constant 
matter energy density, and hence total energy density 
p. The definition of the last scattering hypersurface is 
critical here: note, in particular, that it does not con- 
tain colder and hotter regions which contribute to the 
anisotropics observed at later times, as some accounts 
state (see, e.g., [6, 24, 25]) [26]. Rather, it is a hypersur- 
face of constant (comoving) temperature. Nevertheless, 
in a realistic cosmology, Els is still a perturbed surface, 
with generally non-vanishing intrinsic and extrinsic cur- 
vature and matter perturbations apart from Pm, and 
these perturbations will source anisotropies. Also, be- 
cause of these perturbations, the comoving congruence 
will not in general be orthogonal to Els ■ 

In order to calculate the anisotropies observed at 
late times we must propagate the radiation along null 
geodesies from the observer's LSS. Consider a light ray 
following a null geodesic O with tangent vector that 
extends from a point of emission, E, on the observer's 
LSS to a point of reception, R, and also consider a time- 
like congruence defined in the vicinity of O. We define 
the congruence in order to provide a frame with respect 



to which a local energy density (and hence temperature) 
can be expressed. Then we can decompose i> M at each 
point on O into parts parallel and orthogonal to u M ac- 
cording to 

«" = 7(tt" + n"), (13) 

with n^iifj, = and n M n M = 1. The spatial vector de- 
fines the spatial direction of propagation of the light ray 
at each point along O. Figure 1 illustrates this geometry. 
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FIG. 1: A conformal spacetime diagram showing a light ray 
O emitted from the LSS on Els at event E and received at 
7?. A foliation E t is indicated together with its orthogonal 
congruence tt M , which is comoving on the LSS and at R. For 
clarity, orthogonal vectors are displayed as if the geometry 
was Euclidean. 



Since the radiation is emitted from E with a thermal 
spectrum, an observer at any point on O, with four- 
velocity it M , will observe the radiation travelling along 
O to also have a thermal spectrum with temperature 

T oc -u^Vf, = 7. (14) 

Because the temperature at emission, Te, is defined with 
respect to the plasma frame and is constant on Els, the 
calculation will be simplest if we choose u' 1 to be comov- 
ing with the plasma on the observer's LSS. Similarly, it 
will be most natural to choose it' 1 to be comoving at R, 
since any observer will presumably be composed of mat- 
ter and comoving with it. Note that such an observer, 
comoving with the total (effectively matter) energy den- 
sity at R, will generally not be comoving with the ra- 
diation, and so will observe a dipole anisotropy. As we 
will see, the congruence can be freely chosen in between 
E and R, although some choices will be computationally 
more efficient in practice. In Sec. Ill B 1 we will relax the 
constraint that u M be comoving at E and R. 

Now we will derive an expression for the evolution of 
the "redshift parameter" 7 along O, which will tell us 
how the observed temperature evolves, using only the 
geodesic equation, 

«".„«" = 0. (15) 
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Consider the quantity 



1 



(16) 



= -9lV JlV + (l /. ( 17 ) 

which is related to the expansion rate of the congru- 
ence projected into the plane defined by it^ and n M 
[27] (indeed H,, reduces to the familiar Hubble rate for 
the comoving congruence in a homogeneous and isotropic 
spacetime). Now, using Eqs. (14) and (15), we have 



T 

dX ' 



(18) 
(19) 



where A is an affine parameter along O (note that this last 
expression is independent of the affine parameter chosen) . 
The expression Eq. (19) gives the exact evolution of the 
redshift parameter 7 along O with respect to the con- 
gruence m m in an arbitrary spacetime. It describes the 
familiar behaviour of increasing redshift (decreasing 7) 
in an expanding universe, with "expansion" now seen to 
mean precisely that H n > 0. 

We can in principle integrate Eq. (19) along O to find 
the temperature T^(n M ) observed at R in direction — n M , 
i.e. the CMB temperature sky map. With an appropri- 
ate choice of affine parameter [setting the proportionality 
constant equal to unity in Eq. (14)], we have [28] 









T R (nn = (Tz 1 4 


/ H n dX) . 


(20) 









However, in practice it would be very difficult to deter- 
mine H n as a function of affine parameter along each null 
geodesic in order to perform the integral in Eq. (20) for 
a particular spacetime. Instead, if we define a time coor- 
dinate t by foliating the spacetime into spacelike hyper- 
surfaces of constant t, S t , (see Fig. 1) we can transform 
the integral into one that is more tractable in terms of 
the new coordinate. The most convenient choice for the 
foliation is that which is everywhere orthogonal to the 
congruence u M [29]. 

Using Eq. (13), we can write 



dX 



7 



dT 

N dt 

7. d^' 1 ) 
N dt 



(21) 
(22) 
(23) 



Here r is proper time, N is the lapse function for the 
slicing S t , and the subscript or u M indicates the direc- 
tion in which the derivative is taken. In the intermediate 



steps we have chosen to define 7 away from O to be con- 
stant along the S t . Integrating Eq. (23) along O and 
using Eq. (19) and the proportionality Eq. (14) between 
7 and temperature, we finally obtain 



exp 




H n Ndt 



(24) 



where i£;(n M ) is the value of t at the point of emission E 
on the LSS corresponding to the observed direction — 
at R, and tu is the time of observation. [Note that in gen- 
eral Els w ill n °t coincide with one of the slices Sj; hence 
the dependence t^n^).] This remarkably simple expres- 
sion is exact, and so is not restricted to linear, adiabatic, 
or scalar fluctuations, and it applies to all scales (subject 
of course to our basic assumption of abrupt recombina- 
tion). In particular, Eq. (24) encapsulates in principle the 
acoustic peak structure of the CMB, any ISW contribu- 
tion that may arise, as well as the effect of gravitational 
waves. This equation is a purely geometrical result, in- 
dependent of any dynamical input such as stress-energy 
conservation or Einstein's equations. To our knowledge 
Eq. (24) has not been written down before, although a 
related expression appears in Ref. [7], and a related lin- 
earized expression appears in Ref. [30]. 

Equation (24) tells us that the observed temperature 
in some direction on the sky is determined entirely by the 
integrated line-of-sight component of expansion (or num- 
ber of "e-folds") along the null path from the LSS, of a 
congruence that is comoving with the plasma at the LSS 
and matches the observer's four-velocity at R. Thus we 
can interpret the observed anisotropic CMB sky as a uni- 
form temperature surface viewed through an anisotropi- 
cally expanding universe: the observed hot and cold spots 
on the sky are simply "closer" and "farther" , respectively, 
from us, in terms of e-folds of expansion. Note, however, 
that we cannot view the redshifting as uniquely defined 
at intermediate points between E and R, since we are 
free to deform the congruence between those endpoints. 
Rather, it is only the total integral that is independent of 
the choice of congruence m m between the endpoints, when 
we fix the position and state of motion of the observer at 
R. For the special case of a homogeneous and isotropic 
FRW cosmology, and choosing the comoving congruence, 
for which H n = H and N = 1, we immediately recover 
from Eq. (24) the familiar result for the cosmological red- 
shift 



Tr 
T e 

for scale factor a. 



exp 



Hdt 



a-E 
or' 



(25) 



B. Linearized results 



1. Arbitrary gauge 



The result Eq. (24) is very general, but probably has 
limited direct use for calculating anisotropics. However, 
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it can be straightforwardly expanded to linear (or even 
higher) order in perturbation theory, and such a lin- 
earized calculation will be very convenient, as linear the- 
ory captures very well the evolution of structure at early 
times and very large scales today. In doing this it will 
prove helpful to generalize the result to congruences 
that are non-comoving at E and R. With such thrcad- 
ings it will then be necessary to provide explicit boosts 
at the LSS and at the reception point R to compensate. 
Similarly, at linear order it will be simple to write the in- 
tegral to the LSS in terms of an integral to some constant 
time slice, plus a contribution due to the linear tempo- 
ral displacement to the actual LSS. The boost at the 
LSS will constitute what is often termed the "Doppler" 
or "dipole" contribution to the anisotropies, while the 
temporal displacement contributes to what is sometimes 
called the "monopole" contribution. 

To calculate the boosts, consider the general non- 
comoving congruence and the direction u M comoving 
with the plasma at emission point E (see Fig. 2). If we 
temporarily construct scalar fields t and t in the vicinity 
of E such that = -t^ and ii^ = with t = t-St B , 
then we have 



(26) 



While this expression is exact, the "boost displacement" 
Sts evaluated at linear order is simply the linear tem- 
poral displacement between hypersurfaces orthogonal to 
it M and u^, in the vicinity of E (see Fig. 2); this can be 
readily calculated in the metric formalism as the gauge 
transformation required to take the gauge specified by 
the slicing orthogonal to m m into the plasma-comoving 
gauge. To calculate the change in observed temperature 
due to the boost, we require the quantity 

7=-2%=7(l-V*fl). ( 27 ) 

which is valid at first order. To derive this expression 
we have used the fact that u^St'^ vanishes as first or- 
der, which follows from the normalization of the four- 
velocities. Note that Eq. (27) simply describes a local 
Lorentz transformation of photon energy. This expres- 
sion can also be applied at the reception point R, al- 
though the direction is free in principle there. If we 
choose the observer to be comoving with matter then the 
displacement Sts at R will be given by the gauge trans- 
formation required to take the gauge specified by into 
the comoving gauge. The freedom to choose at R only 
effects the dipole anisotropy at linear order, as shown in 
Appendix C. 

To calculate the temporal displacements at E and R, 
we can write the integral in the exact expression Eq. (24) 
as 



H n Ndt= / H n Ndt, (28) 

t R Jt R +8t D (R) 



where and label particular slices E t , and can be 
considered the background emission and reception times. 
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FIG. 2: A conformal spacetime diagram showing a light ray 
O emitted from the LSS at E. The arbitrary foliation Ej is 
indicated together with its orthogonal congruence . Also, 
the hypersurface E 9 orthogonal to the direction {t M comov- 
ing with the plasma is indicated, together with the "boost 
displacement" StB, which takes E t „(„M) into E 9 . 



The displacement St o (E) accounts for the separation be- 
tween the background slice Ej E and the true last scatter- 
ing hypersurface Els, and is a function of n M . Similarly, 
the displacement St p (R) accounts for the separation be- 
tween the background slice Ej„ and the actual slice on 
which the reception point R is located. Sto{R) can be 
considered a function of position if we wish to evaluate 
the anisotropies at various reception points R. At linear 
order, we can then write 



t E (n") 



H n Ndt 



H n Ndt+ (HSt D ) 



(29) 



where H (t) is the zeroth order (background) Hubble rate 
and we have assumed that the lapse equals unity at ze- 
roth order, so that the coordinate t is a perturbed proper 
time. 

The displacement St^iE), like the boost displace- 
ments, can be readily calculated in the metric formal- 
ism as the gauge transformation required to take the 
gauge specified by the slices E t (orthogonal to u M ) into 
the gauge specified by Els, i- e - the uniform radiation en- 
ergy density gauge. The displacement Sto(R) is not fixed 
uniquely by any such physical prescription, but only af- 
fects the monopole anisotropy, as shown in Appendix C. 

The temperature observed at R in direction —n^ can 
now be written using the boost relation Eq. (27) as 



. 1 - njt'fi , (30) 

at linear order. The prefactor on the right-hand side of 
this equation is simply the redshift due to the expan- 
sion along the congruence u M , so it is given by Eq. (24). 
Writing 



H n N = H + S(H n N) 



(31) 
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and using Eq. (29), we have 



— - — - = cxp 



rt E (n») 



H„Ndt 



(32) 



^ (l + JJ 5(H n N)dt+ (H5t D )\ E \(33) 



at linear order, where we have defined the "background" 
observed temperature Tr by 



(34) 



Finally, combining Eqs. (33) and (30), and defining 
5T(n») = T R (n») - f R , we can write 




Tr 



S(H n N)dt + (H6t D 



(35) 



Equation (35) gives the observed temperature 
anisotropy at linear order in terms of the linc-of-sight 
expansion perturbation S(H n N) in an arbitrary gauge 
(specified by the hypersurfaces orthogonal to u M ) and the 
temporal displacements Sts and 8to required to trans- 
form from that arbitrary gauge to comoving and uni- 
form density gauges. The only approximations involved 
in Eq. (35) are those of abrupt recombination and lin- 
earization. The terms HSto and n^Stg evaluated at the 
LSS are sometimes called the "monopole" and "Dopplcr" 
contributions, respectively. The geometrical nature of 
the terms in Eq. (35) provides a clear and unambigu- 
ous interpretation of the anisotropy, without reliance on 
coordinate-dependent notions such as gravitational po- 
tentials (see Rcf. [31] for a related discussion). 

Equation (35) is in a form that makes it easy to evalu- 
ate the anisotropy using any gauge for the perturbations 
that we choose. First, for the line-of-sight integral, by 
linearizing the exact expression Eq. (17) we have 



6(H n N) = -89 + a^n^n" + + HSN. (36) 
o 

The geometrical quantities in this expression can be writ- 
ten in terms of the metric perturbations using Eqs. (A2), 
(A5), (A6), and (A8), giving 

S(H n N) = -j, + ^n* 1 + a, liV n' 1 n v + ifZ^n'W. (37) 



P(y), bpMi an d 9(7) i while at the observation point 
R the transformations are determined by the hypersur- 
face and state of motion of the observer chosen, as we 
will see. 

Since the quantity T/?(n M ) is observable, the general 
expression Eq. (35) must be independent of the gauge or 
congruence chosen, if the point R and the four-velocity 
of the observer are held constant. This is demonstrated 
explicitly in Appendix B. However, if the observation 
point and four-velocity are allowed to transform with the 
gauge transformation, then the anisotropies will depend 
on the gauge, as shown in Appendix C. Expanding the 
anisotropy in terms of the spherical harmonics, Y^ m (n /i ), 
the multipole amplitudes are 



Ofn 



ST(n^ 
Tr 



-y/ m K)dfi. 



(38) 



We show explicitly in Appendix C that only the dipole 
anisotropy a\ m and monopole perturbation aoo change in 
the latter case, at linear order. 



2. Recovering previous results in zero-shear or longitudinal 

gauge 

We will now illustrate the usefulness of the general 
linear expression for the anisotropies, Eq. (35), by calcu- 
lating in a particular gauge the anisotropy due to both 
adiabatic scalar and tensor sources, recovering previous 
results. We require both the line-of-sight integral in that 
expression as well as the temporal displacements St jj and 
5ts for the boundary terms. If we choose the congruence 
such that the scalar-derived part of the shear van- 
ishes at linear order, then the integrand Eq. (37) takes a 
particularly simple form. The frequently used longitudi- 
nal gauge has this property, giving 

S(H n N)^-ij a + ^n^ + ^H^n^. (39) 

In these expressions, the subscript a indicates a zero- 
shear or longitudinal gauge quantity, and the overdot in- 
dicates the proper time derivative in the direction of . 
Therefore, using the first order expression 



dt 



(40) 



we can write the integral in Eq. (35) as 



Here ^ is the curvature perturbation, <fi is the lapse per- 
turbation, a is the shear scalar, and H^ v is the tensor 
metric perturbation. Equation (37) is valid in arbitrary 
gauges, i.e. for arbitrary congruences m' 1 , and it is now 
trivial to fix the gauge, as we will see in Sec. Ill B 2. Sec- 
ond, for the boundary terms in Eq. (35), we can work 
out the required gauge transformations Stu and Sts us- 
ing Eq. (A14) and (A15). At the emission point E those 
transformations arc applied to the radiation quantities 



S(H n N)dt 



4> 



1p(7 + <Pa 



-H^n v \ dt 



(41) 



It is simple to calculate the displacements Stu and Sts 
for the case of large-scale adiabatic modes, for which the 
uniform radiation energy density and uniform total en- 
ergy density hypersurfaces coincide, and for which the 
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plasma-comoving and total comoving directions coincide. 
(Additionally, these surfaces and directions coincide even 
on small scales for the artificial case of a cold dark mat- 
ter free universe, in the tight coupling approximation.) 
Although Els is defined as a surface of constant pM, it 
will be easier to calculate the position of uniform total 
density surfaces. The temporal displacement that takes 
us from zero-shear to uniform total energy density gauge 
is, using Eq. (A14) and the linearized energy constraint 
equation Eq. (A16) for the total energy density pertur- 
bation 5p, 



5t D = - 



Sp a 



3H U a + H<t> a ) - ivV 



12irGH(p + P) 



(42) 



This expression can be applied at point E on the LSS as 
well as at the reception point R if we choose to place R 
on a surface of uniform energy density. The "boost dis- 
placement" that transforms from zero-shear to comoving 
gauge is, using Eq. (A15) and the linearized momentum 
constraint equation Eq. (A17) for the total momentum 
density scalar q, 



StR = - 



Ipcr + H4>a 

4nG(p + P) ' 



(43) 



Again, this can be applied at both E and at R if we 
choose the observer to be comoving. 

Equations (41) to (43) with Eq. (35) completely specify 
the anisotropics in terms of quantities in zero-shear or 
longitudinal gauge. In practice it is common to make 
approximations. If we assume that the anisotropic stress 
is negligible (as is the case for matter or A domination), 
then we have = </v (see, e.g., [32]). As explained 
above, the displacements 8to and Sts at the observation 
point R only affect the observed monopole and dipolc. 
Dropping these terms, and using the background energy 
constraint, the anisotropy for I > 1 then becomes 



tB 



Wo 



2 H 2 



2 1 
3#2 



(44) 



All quantities outside the integral are to be evaluated at 
the point on the LSS corresponding to viewing direction 
—n^. This expression agrees precisely with Eq. (4.7) in 
Ref. [8], for the case of adiabatic perturbations, including 
even a term the authors of [8] describe as arising from 
"subtle gauge effects." 

Making further simplifications, we have ip a = at 
last scattering if the pressure vanishes exactly there. If 
we consider anisotropies on the largest angular scales, 
sourced by modes with comoving wavenumber k <C aH , 
we can drop all the gradient terms in Eq. (44). The result 



5T(n^) 



Wo 



-H^n" ) dt+-ME), (45) 



in agreement with the well-known result for the Sachs- 
Wolfe effect due to large-scale scalar and tensor sources. 
Note that this result includes a part that is evaluated at 
the boundary E, which comes from both the temporal 
displacement Stu(E) and the integral in Eq. (35). The 
remainder of that integral, which cannot be placed at 
the boundary, appears in Eq. (45) as a contribution that 
is due to physical metric fluctuations along the line of 
sight. The scalar part of this contribution is known as the 
integrated Sachs- Wolfe effect. During matter domination 
tp a = 0, so all scalar effects of the perturbed expansion 
along the line of sight can be placed at the boundary. 



IV. TRANSFER FUNCTIONS 
A. Dipole 

Next we will use the results derived so far to per- 
form a careful calculation of the CMB dipole due to 
scalar sources. More precisely we will calculate the dipole 
power, or the variance in the dipole anisotropy ai m , 
namely 



Ci 



l(ai m )| 2 , 



(46) 



over realizations of the assumed Gaussian random pri- 
mordial fluctuations. (The independence of C\ on m will 
follow from statistical isotropy.) For the dipole the choice 
of frame for the observer at R is critical, since a boost 
at the observation point changes the dipole according to 
Eq. (C12). For example, if we choose the observer's frame 
to be comoving with radiation, so that the radiation mo- 
mentum density q^ vanishes, then trivially the observed 
dipole vanishes. Indeed, combining Eqs. (A12) and (C12) 
we have 



4p 2 



(7) 



(47) 



so that the magnitude of the observed dipole in any 
frame is proportional to the radiation flux observed in 
that frame. We will adopt the most natural and phys- 
ically best-motivated choice, namely the frame comov- 
ing with matter (essentially the total comoving frame at 
late times) for the calculation of C\. With this choice 
of frame, Eq. (47) has the simple interpretation that the 
observed dipole is a measure of how well worldlincs co- 
moving with radiation and matter coincide. 

We begin with the general linear result, Eq. (35). All 
of the parts of this equation were carefully calculated in 
zero-shear gauge for a comoving observer in Sec. IIIB 2. 
We can ignore the monopole contribution St u (R) since it 
only affects the I — mode. The displacements 6tr>(E) 
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and 8tB were calculated in Eqs. (42) and (43). Com- tcring (so ipa(E) = 0), and finally ignoring the term 

billing these results with Eq. (41) for the line-of-sight -t-V 2 ^v(.E), we have 
integral, and ignoring anisotropic stress at all times (so 
i\) a = assuming matter domination at last scat- 



«v> _ \ ME) _ pM^i + _ 1 \ ^mn + 2 f ^, (48) 



Tr 3 3 \3' / H R j E 

Here we have used Eq. (A21) for the growth function g R = g(tu) and the relation Eq. (A24) to simplify the expression. 
The approximation of matter domination at last scattering (which implies zero anisotropic stress) results in errors in 
Gi for small i on the order of 10% [33]. Neglecting the term ^ 7 V 2, (T (i?) is entirely justified considering that Eq. (42) 
for St u (E) employed the approximation that the uniform radiation energy density and uniform total energy density 
hypersurfaces coincide, which is only valid on large scales (scales that were super- Hubble at last scattering). The 
contribution to the dipole from last scattering will be dominated by these large scales, as wc will see. 

To calculate the variance of the dipole, it will be helpful to expand the function tjj (T (E) in spherical harmonics as 



/ "ftwW^RL^N^K)' ( 49 ) 



Here ji is the spherical Bessel function of the first kind, tls is the comoving radius of the LSS, k is the comoving 
wave number, and T(k) is the transfer function defined in Eq. (A19). With the aim of expressing C\ in terms of 
the primordial spectrum Vn(k), we have used Eq. (A22) to write i\) a in terms of the primordial comoving curvature 
perturbation lZ pr . Similarly we will need the expansion of the "Doppler" contributions at E and R, 

n^(E) = -L J dkk 2 T(k)Y,KL( k )ji( k ^s)Y em (~nn, (50) 

and 

n^\R) = lim iJl^ [ dkk 2 T(k) Y>£„(A0iK (51) 

= \J^— I dkk 2 T(k)J2K r m (k)Y lm (-nn, (52) 
5 V 7r an J 



where the prime indicates differentiation with respect to the argument, and we have used the relation ^(0) = 5}/3. 
Inserting these expressions into Eq. (48), we can evaluate the dipole anisotropy using Eq. (38) (with £ = 1) and the 
orthonormality of the Y^ m . This gives 

a lm = -\J- I ' dkkKf^T^Txik), (53) 



5 V 7T _ 

where T\{k) is a new transfer function, called the dipole transfer function, and defined by 



Ik ( 5\ k f R 

Ti(fc) ^ jiCfcr-Ls) + _- 7 ;(fcr LS ) +(g R --z) —7T + 6 / g(t)h[kr(t)]dt. (54) 

dEtlE \ oj OrUr Je 



Here r(t) is the comoving radial coordinate of the light 
ray as a function of the time coordinate. Finally, using 
Eq. (46) and the statistical relation Eq. (A25) for the 
primordial power spectrum Vn(k), we find for the dipole 



power 

Cl = 1/ ^n(k)THk)Tl{k). (55) 

Note the third term in Eq. (54), which is proportional 
to k and dominates on small scales today. It might ap- 
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pear that this term would lead to an ultraviolet diver- 
gence for the dipole, but Eq. (55) is rendered finite by 
the function T(k), which decays like k~ 2 on small scales. 
This results in a dipole amplitude \[C\ ~ 10 -3 for stan- 
dard cosmological parameters. 

B. Monopole 

The monopole temperature perturbation is the I = 
component of Eq. (38), i.e. 

(56) 

and its variance over realizations of the primordial fluctu- 
ations will be called Co- As discussed above Eq. (B14), 
the monopole so defined is ambiguous in that different 
choices of the background times tg and in lead to dif- 
ferent etoo- This well-known freedom is equivalent to our 
inability to uniquely separate a "background" CMB tem- 
perature from the observed mean temperature, which 
would be required to define a monopole perturbation. 
Hence the monopole perturbation in Eq. (56) cannot be 
observable. 

Nevertheless, it is still possible to sensibly define a 
monopole and consider its theoretical properties. Recall 
that for the dipole it was necessary to fix the observer's 
frame by a physical prescription (namely that it be co- 
moving with matter) . Analogously, to define a monopole 
the key point is that we must fix, by a physical prescrip- 
tion, the spacelike hypersurface on which we place the 
observer. Clearly there is freedom in how wc choose this 
slice. Recall that an observer chosen to comove with ra- 
diation observes no dipole. The analogous situation with 
the monopole is an observer placed on a uniform radia- 
tion energy density hypersurface, for which the calculated 
Co must vanish. Analogously to Eq. (47) for the dipole, 
we can write 

l«oo| 2 = j^H^( 7) ) 2 , (57) 

4/? (7) 

where 5pM is the radiation energy density perturbation 
on the same slicing used to define the monopole aoo- 



Therefore the monopole defined with respect to any slic- 
ing is proportional to the radiation density perturbation 
for the same slicing. 

The simplest and most natural choice of slice on which 
to place the observer is that of uniform matter (essen- 
tially total) energy density. This choice is simplest be- 
cause it requires knowledge of just the local density, p, 
which acts as a clock. It is natural because, as we will 
see, it exhibits close analogy with the dipole. Comov- 
ing slices could also be used, although their construction 
as hypersurfaces orthogonal to the comoving worldlines 
is more elaborate [34]. Thus, while the dipole defined 
above is a measure of the degree to which comoving radi- 
ation and matter worldlines coincide, the monopole de- 
fined here has the simple interpretation as a measure of 
how well hypersurfaces of uniform radiation and mat- 
ter energy density coincide [35]. The arbitrariness in aoo 
mentioned above (due to the freedom to choose the back- 
ground times) only amounts to a constant shift to aoo, 
so its variance is unchanged. This is why, although the 
monopole perturbation aoo is ambiguous and unobserv- 
able for any single observer, the variance Co is still well 
defined theoretically (and could, in principle, be approx- 
imated through observations [36]). 

It should be clear immediately that there is a seri- 
ous problem with attempting to define the variance of 
the CMB temperature on a uniform matter density slice. 
Namely, matter has of course entered the nonlinear stage 
on small scales, and hence hypersurfaces of constant mat- 
ter density cannot actually be defined. Nevertheless, 
smoothing over small scales can recover meaningful lin- 
ear results, at the expense of further ambiguity in the 
form of the smoothing scale. This issue arises because 
in calculating the monopole as defined here, we will see 
that the dominant contribution will be given by the vari- 
ance ((Sp/p) 2 ) evaluated at R, which is dominated by 
small scales (where the slicing is irrelevant). The im- 
portant practical aspect of this brief examination of the 
monopole will actually be in understanding the effects of 
long- wavelength sources, for which the difficulties with 
small scales do not arise. 

Proceeding as with the dipole above, the relevant con- 
tributions to the general linear result Eq. (35) give 



- \ M E) - p*m + (k> - 2) MR) - + 2 f'i>,M. (58) 

3 3 H E \3 n J 90 m a R H R J E 



We can ignore the "Dopplcr" contribution Sts at R since it only affects the t—1 mode, and the same approximations 
have been applied here as for the dipole Eq. (48). 

To calculate the variance of the monopole we will need, in addition to Eqs. (49) and (50), the "monopole" term 
evaluated at the reception point R, 



MR) = -\\-9R f dkkT{k)Y^nf m (k) 3t {Q)Y lm {^) (59) 
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f oo, 



(60) 



= -\^rJ dkkT(k)K%(k)Yt 
using the relation je(0) = 6®. Proceeding as for the dipole case, we find 

fl oo = -\\[\ j dkknZ(k)T(k)To(k), (61) 
where Xq(/c) is the monopole transfer function defined by 

T (k) = j (kr hs ) + —^f (kr hs ) + 5-6g R + -^-—^ + 6 / g(t)j [kr{t)]dt. (62) 



Finally, using the statistical relation Eq. (A25) to evalu- 
ate the variance wc find 



C 



-in 

25 



dk 



V n (k)T 2 (k)T 2 {k). 



(63) 



Note the presence of the Laplacian term (oc k 2 ) in To(fc), 
which will dominate on small scales, and implies that by 
calculating Co we are essentially calculating the (effec- 
tively gauge independent) matter variance ((Sp/p) 2 ) at 
the observation point. In this monopole case the result- 
ing divergence is too strong to be saved by the transfer 
function T(k). Therefore wc expect the variance Co to di- 
verge on small scales, which is simply a reflection of the 
nonlinear nature of matter fluctuations on small scales 
today, as predicted above. That is, the monopole vari- 
ance as we have defined it cannot be quantified. Again, 
the importance of Eq. (62) will lie in its long- wavelength 
behaviour. 



V. LONG- WAVELENGTH BEHAVIOUR 

Recall that all of our approximations have been good 
on very large scales. In this section we examine the long- 
wavelength limit of the Ti(fc) and T (k) transfer func- 
tions. The dipole transfer function for a comoving ob- 
server, Eq. (54), is plotted in Fig. 3, together with the 
monopole function for an observer on a uniform energy 
density slice, Eq. (62), and the transfer functions for 
I = 2, 3, and 4. The transfer functions for I > 1 can 
be calculated from Eq. (48) in exactly the same way as 
for the dipole, with the result 



Te(k) = ji{kr- LS ) + 



2k f R 
-S-J/(*as) + 6 / g{t)ji[kr{t)]dt. 
aE^E Je 

(64) 

Note that the large-scale approximations involved in 
Eq. (48) imply that this expression is only valid for scales 
that are super-Hubble at last scattering. [The transfer 
functions Ti(k) and T (k) are valid for small scales, since 
for large kr^s the second-to-last terms in both Eqs. (54) 
and (62) dominate. These terms are generated locally at 
the observation point R.] A value of flA = 0.77 today 
was used for all of these calculations. 




FIG. 3: Dipole transfer function Ti(k)T(k) for a comoving 
observer (solid curve). For comparison, the monopole transfer 
functions for an observer on a uniform energy density slice, 
To(k)T(k), and the transfer functions for £ = 2, 3, and 4 are 
also shown. Absolute values are plotted. The scales fcrxs = 
10 2 and Atls = 1 correspond roughly to the Hubble scales at 
last scattering and today, respectively. 



Using asymptotic forms for the Bessel functions [37], 
we can show from Eq. (64) that Tg{k) should decay like 
k e as fcrLS ~~ > 0: f° r ^ > 1- This is verified in Fig. 3. 
However, the figure also shows that T\(k) does not decay 
like k for small k; instead it decays like fc 3 , which is faster 
than the decay rate of T2(k). 

To examine the behaviour of the transfer function 
T\(k) in the limit k — > 0, we can use the small-argument 
approximations to the Bessel functions [37] to give 



Ti(k) 



1 



9R ~3 
Oikr^s)' 



OirHr 



3 V " 



g{r))drj 



(65) 



where the conformal time 77, defined via dr] = dt/a, has 
been used to simplify the expression. In writing Eq. (65), 
we have used the approximation g(i_e) = 1, although for 
our numerical calculations we have evaluated Eq. (54) 
without approximation. It is not obvious from Eq. (65) 
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that the O(k) term, in square brackets, vanishes. How- 
ever, the numerical calculation of Fig. 3 shows that the 
0(k) term does indeed vanish, as pointed out above. 

This is illustrated in greater detail in Fig. 4. There we 
have plotted separately the first two terms in Eq. (54), 
which originate at the LSS (called SW^ in the plot), the 
third term in Eq. (54), which comes from the observa- 
tion point R (SWr), and the line-of-sight ISW term. We 
can see that each component separately does scale like k 
on large scales, in particular the LSS term SW^ scales 
like the predicted k , but the sum demonstrates that the 
0(k) terms all cancel. It is the local contribution, SW^, 
not present for £ > 1, which enables the cancellation. 
(Of course the individual SW and ISW components are 
not separately observable!) It is possible to show an- 
alytically this 0(k 3 ) dependence in the special case of 
a cosmological-constant free Einstein-de Sitter universe. 
Then Eq. (65) becomes 



Ti(k) 



30 



(^s+3^ s ^s)+0(fc 5 ). (66) 




FIG. 4: Dipole transfer function Ti(k)T(k) for a comoving 
observer. Absolute values of individual contributions from 
Sachs- Wolfe terms evaluated at emission, E, and at the obser- 
vation point, R, as well as the line-of-sight ISW contribution, 
are indicated. 



Figure 4 also illustrates that the source of the small- 
scale increase in T\{k) is due to the local contribution 
at R, as explained above. The dipole transfer function, 
Eq. (54), suggests a divergence in C\ like a power of k 
on small scales, although this is moderated by the trans- 
fer function T(k). The result is that the observed dipole 
amplitude is a factor ~ lfj 2 larger than the other multi- 
poles. Figure 3 actually illustrates this directly: Eq. (55) 
generalizes to 



Ct = % I T rKik)T2(k)Tm 



(67) 



for all £ > 0. Therefore the expected multipolc amplitude 
is simply proportional to the area under the appropriate 



curve in Fig. 3 [assuming a nearly scale-invariant pri- 
mordial spectrum Viz(k)]. Geometrically, the comoving 
matter and radiation worldlines coincide on large scales, 
and begin to diverge strongly on small scales. The contri- 
bution to the dipole from last scattering, SW#, will vary 
smoothly if the observer's position is varied, whereas the 
local contribution, SW.R, will vary greatly. 

These results indicate that the dipole defined with re- 
spect to the comoving frame receives strongly suppressed 
contributions from super-Hubblc modes. This was in fact 
noted some time ago [15] in relation to the Grishchuk- 
Zel'dovich effect [14], but in the context of a matter- 
dominated universe. In fact, it was claimed in [15] that 
this cancellation would not persist in the presence of 
a cosmological constant, so that super-Hubblc modes 
would have an observable imprint on the CMB, a claim 
which we have now demonstrated to be incorrect. 

As we did for the dipole, we can examine the behaviour 
of the monopole transfer function To(k) in the limit k — » 
0, giving 



T (k) 



1 



'LS 



2r LS 
So-eHe 



2<7fl 



3a 2 E H E 



g{t)r 2 {t)dt 



0{k A 



(68) 



where fi m is the standard matter density parameter. In 
this monopole case, we have thus shown analytically that 
the leading-order terms in the transfer function. Eq. (62), 
which go like k°, do cancel in the presence of a cosmo- 
logical constant. However, Fig. 3 shows that an even 
stronger result holds in the monopole case: the 0(k 2 ) 
terms in Eq. (68) actually cancel as well! Again, this is 
shown in greater detail in Fig. 5, where it is apparent 
that an exquisite cancellation occurs in the individual 
components, which scale like k° on large scales, to give 
the 0(k 4 ) total transfer function on large scales. Fig- 
ure 5 also illustrates the local source, SW,r, of the strong 
small-scale divergence discussed above. Geometrically 
the hypcrsurfaces of uniform matter and radiation den- 
sity coincide on large scales, and begin to very strongly 
diverge on small scales. This divergence makes it impos- 
sible to quantify the total power C'o with our definition 
of the monopole. 



VI. DISCUSSION 

Our result in Sec. V that the dipole and monopole 
receive suppressed contributions from large scales, even 
in the presence of a dominant cosmological constant, 
strongly suggests that the cancellations involved are not 
accidental. To understand the origin of this behaviour, 
consider first the monopole case. Physically, the suppres- 
sion as k — > in the monopole transfer function shown 
in Figs. 3 and 5 means that surfaces of uniform total and 
radiation energy density coincide on the largest scales, 
according to our definition of the monopole in Sec. IV B. 
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FIG. 5: Monopole transfer function To(k)T(k), for an ob- 
server on a uniform energy density slice. Absolute values of 
individual contributions from Sachs- Wolfe terms evaluated at 
emission, E, and at the observation point, R, as well as the 
line-of-sight IS W contribution, are indicated. 



But this is just the statement of adiabaticity: an adia- 
batic matter-radiation fluid is characterized by the con- 
dition 



Sp 
P 



(69) 



on the largest scales. Therefore, Eq. (A14) shows that the 
same gauge transformation takes us to both constant to- 
tal matter and constant radiation hypersurfaces; in other 
words, those surfaces must coincide. It is important to 
point out that adiabatic long-wavelength modes remain 
adiabatic under evolution [32], and indeed the comoving 
curvature perturbation 1Z remains constant in time (see, 
e.g., [38]). This trivial evolution on supcr-Hubblc scales 
means that when constant matter and radiation density 
surfaces coincide on large scales at last scattering, due 
to adiabatic initial conditions, they must also coincide 
today. 

An analogous situation holds for the dipolc. In this 
case, the adiabaticity condition implies that, on large 
scales, the radiation and total matter comoving world- 
lines coincide (i.e., there is no "peculiar velocity" isocur- 
vature mode between the two components). According 
to our definition of the dipole in Sec. IV A, this is sim- 
ply the statement that the dipole is suppressed on large 
scales. 

This leads us to the important conclusion that this 
insensitivity to long-wavelength sources must apply re- 
gardless of the matter content of the Universe, as long as 
adiabaticity holds. Therefore the suppression we found 
for the specific case of cosmological constant (or Einstein- 
de Sitter) universes must in fact occur in general. Note 
that this may have relevance to very recent discussions 
regarding a potential power asymmetry in the CMB 
[39]. The exquisite cancellations visible at large scales 



in Figs. 4 and 5 between the SW and ISW components 
illustrate a previously unrecognized relation between the 
two, which is enforced by the condition of adiabaticity. 

In brief, the dipole and monopole as defined here are 
just measures of departures from the adiabaticity con- 
dition, Eq. (69), which generally occur on small scales. 
Are these results sensitive to the definitions of dipole and 
monopole used? As long as the dipole and monopole are 
defined physically, i.e. in relation to locally mcasureable 
quantities, then the results must still hold. An impor- 
tant example which does not satisfy this criterion is the 
zero-shear, or longitudinal gauge, since it is not defined 
in terms of local, observable matter quantities. If the 
dipole (or monopole) is defined with respect to a zero- 
shear frame, then it may appear from theoretical cal- 
culations that the dipole is sensitive to long-wavelength 
modes. However, this sensitivity cannot be observable, 
since zero-shear frames cannot be uniquely constructed 
locally. (A linear boost or "tilt" of a zero-shear frame is 
still a zero-shear frame.) Thus great care must be taken 
when considering behaviour on large scales with longitu- 
dinal gauge. 

Indeed, another way to understand this result is to re- 
alize that, in the limit k/(aH) — ► 0, an adiabatic pertur- 
bation mode locally becomes essentially pure gauge, and 
can be removed within any sub-Hubble region by a sim- 
ple boost coordinate transformation. Such a mode is in- 
distinguishable locally from a homogeneous background, 
and hence cannot have any observational consequences 
such as a dipole anisotropy [40]. 

Finally, we note that when the condition of adiabatic- 
ity is relaxed, then our conclusions no longer hold [16]. 
In the presence of isocurvature perturbations, it is possi- 
ble that a physically defined dipole be sensitive to super- 
Hubble modes, since the extra freedom allows for a rela- 
tive tilt between comoving matter and radiation comov- 
ing worldlines on large scales. 

While we have emphasized here the consequences of 
adiabaticity, it is hoped that other applications will fol- 
low from the exact formalism for CMB anisotropics that 
we have developed. One possibility is the evaluation of 
anisotropics, in particular the ISW effect, in void mod- 
els of acceleration [41] (see [42] for a brief review), which 
have not yet been confronted with observations at the 
perturbative level [43]. Note however that the present 
approach is not limited to calculating CMB anisotropies, 
and that more generally it is applicable to calculating 
redshifts in arbitrary spacetimcs. Potential uses include 
calculating the redshift-luminosity distance relation in 
perturbed spacetimes (see, e.g., [44]). 

Note added: When this work was essentially complete, 
a related paper appeared [45] , which appears to support 
our conclusion that long-wavelength perturbations can- 
not effect the CMB dipole. 
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APPENDIX A: METRIC-BASED APPROACH 

In this appendix, we will collect together several ele- 
ments of the metric-based approach to cosmological per- 
turbations that will be useful in describing the CMB 
anisotropies at linear order (see, e.g., [32] for a review 
of metric-based perturbation theory) . In the metric ap- 
proach a set of coordinates are defined in the spacetime 
by a foliation into spacelikc hypcrsurfaces of constant t, 
£t, and a threading into timelike worldlines of constant 
spatial coordinates x 1 , where latin indices run from 1 to 
3. The gauge freedom of perturbation theory is related to 
our freedom to choose such a slicing and threading. How- 
ever, it can be shown (see [46, 47]) that at linear order, 
physical perturbations in FRW backgrounds are gauge 
invariant under changes in the threading. Therefore the 
gravitational dynamics can be entirely expressed in terms 
of spatially gauge-invariant quantities. The reason for 
this invariance is simply the homogeneity of the back- 
ground spacetime, and its importance is that it means 
that there is effectively just a single degree of gauge free- 
dom on FRW backgrounds, namely the temporal position 
of St at each event. Thus, to simplify expressions to fol- 
low, we will choose the congruence of coordinate threads 
(with tangents u^) to be orthogonal to the S t , so that the 
shift vector (metric component goi) vanishes. These spa- 
tial coordinates are comoving at zeroth order, but may 
depart from comoving at first order. With this choice, 
there is a direct correspondence between the gauge free- 
dom (i.e. the freedom to choose the slicing) in the met- 
ric formalism, and the freedom to choose the congruence 
in the covariant formalism. Evaluating the metric using 
Eq. (1) in the chosen coordinates we have 

,9oo = -N 2 , g lj ^h ij , (Al) 

where N is the lapse function for the slicing. We con- 
sider scalar and tensor perturbations only, as vectors are 
ordinarily thought to be cosmologically irrelevant. 

The spacetime is completely described by the lapse 
function, the intrinsic curvature of the spatial metric h^, 
and the extrinsic curvature of the S t . We define the lapse 
perturbation <f> through 

N = 1 + 4>. (A2) 

The only part of the intrinsic curvature of h^ v that we 
will need is the perturbed Ricci scalar, S^ 3 'R, for the spa- 
tial slices S t , which defines the curvature perturbation ip 
via the relation 

S^R = 4rVV (A3) 



Here V 2 /a 2 = D^D^, for background scale factor a, so 
V 2 is the comoving Laplacian. The extrinsic curvature 
of the slicing is specified by the expansion and shear of 
the normal congruence to the S t , and is related to the 
spatial metric through [22] 

^9h ij +a ij = ^h ij , (A4) 

where the overdot represents the proper time derivative 
along u M . At linear order the trace of this equation gives 
for the expansion perturbation 

56 = -3H(t> -3ip + -^V 2 cr, (A5) 

where H is the background Hubble rate. The shear scalar 
a describes the scalar-derived part of the shear, with 

(Tij = D.Dja - ^2^ 2(Th v + ( A6 ) 

where is the transverse and traceless (tensor) part of 
hij , and we have ignored the vector-derived part of the 
shear. One further important quantity is the acceleration 
of the worldlines normal to the S t , Eq. (4). It is related 
to the lapse through the exact expression 

a„ = ^D^N, (A7) 
which at linear order becomes 

a M = D^. (A8) 

The various quantities defined above can be related to 
the explicit component form of the metric at linear order 
(ignoring the vector part), 

Sg 00 = -2(f>, S gij = a 2 (-2^7« + 2S,y + Hij), (A9) 

where jij is the background spatial metric, and the trace- 
free scalar part of Eq. (A4) gives a = a 2 E. 

Once the slicing S t is specified, i.e. the time coordinate 
is chosen, then the perturbation in any exact quantity 
X(x l , t) is fixed via 

5X{x\t) = X(x\t) - X(t), (A10) 

where X(t) is the homogeneous background value. 
Therefore our freedom to vary the S t (equivalently to 
vary the orthogonal congruence u^) results in an inher- 
ent ambiguity in our ability to specify any perturbation. 
This temporal gauge freedom can be used to simplify cal- 
culations by choosing an appropriate congruence, as we 
see with the Sachs- Wolfe calculation in Sec. IIIB2. It is 
straightforward to calculate the change in any perturba- 
tion under a gauge transformation t — > t — St (see, e.g., 
[47]). A few results that we will need are, at linear order, 

5p^Sp + pSt, (All) 
q^q-(p + P)St, (A12) 
tp -»■ V> - HSt, (A13) 
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where q is the scalar from which the momentum density 
q^ L is derived, through = D^q. These results imply 
that the gauge transformation required to go from an 
arbitrary initial gauge to uniform energy density gauge, 
defined by Sp = 0, is 



st = - s 4, 
p 



(A14) 



and the transformation that takes an arbitrary initial 
gauge into comoving gauge, defined by q = 0, is 



St 



'1 



P 



(A15) 



Unless otherwise stated, all expressions in this work will 
be presented in an unspecified gauge, i.e. they will apply 
to arbitrary gauges. 

We will also need the Einstein constraint equations in 
order to relate the matter perturbations to the metric 
perturbations. Projecting the Einstein equations twice 
along and linearizing, we find the energy constraint, 

m(ip + H<p) - -jV 2 (ib + Ha) = -AirGSp. (A16) 

Similarly, projecting once with and once with h^ v 
gives the linearized (scalar) momentum constraint 



ip + H(f> = -AirGq. 



(A17) 



It is conventional to express the primordial power spec- 
trum in terms of the comoving gauge curvature pertur- 
bation, usually denoted 1Z, since 7Z is constant on large 
(super-Hubblc) scales for adiabatic modes and hence its 
value at late times can be trivially related to the predic- 
tions of an inflationary model (see, e.g., [38]). However, 
we will see that it is simplest to perform the linear Sachs- 
Wolfe calculation in terms of the zero-shear gauge curva- 
ture perturbation, ip a , where we denote zero-shear gauge 
perturbations by the subscript a. For matter domina- 
tion, it is simple to relate ip a to 1Z by performing a gauge 
transformation from zero-shear to comoving gauge; using 
Eqs. (A12) and (A13) gives 



:1Z. 



(A18) 



by the presence of its argument, k.) The transfer function 
T(k) approaches unity at small fc, and decays roughly like 
k~ 2 for k > fc eq , where k eq is the wavenumber which en- 
ters the Hubble radius at matter-radiation equality (see, 
e.g., [38]). T(k) can be calculated accurately numerically, 
e.g. using packages such as CAMB [10, 11]. Combining 
Eqs. (A18) and (A19) gives 



--T(fc)ft pr (k). 
5 



(A20) 



The perturbation is constant during matter dom- 
ination, but it decays once a cosmological constant be- 
comes important. This decay is independent of scale and 
can be described by a function g{t) via 



which gives 



ip a {k,t) 



--g(t)T(k)K^(k). 



(A21) 



(A22) 



The function git) approaches unity and zero at early and 
late times, respectively, and solving the linearized dy- 
namical Einstein equation for ip a gives 



5 Q, m oH H 

g(t) = ^ — - — 

2 a 



dt' 



(A23) 



Here Hq is the background Hubble rate today, and £l m .o 
is the ratio of matter to total energy density today. Using 
this last expression we can derive the useful relation 



— + q = s2 r 
H 



2 2 9 



(A24) 



The statistics of the assumed Gaussian random primor- 
dial fluctuations are completely described by the power 
spectrum Vn(k), defined by 

(ft pr (k)ftP"(k')) = 2ir 2 5 3 (k - k')^P. (A25) 

k 

With this definition Vn(k) is dimensionless, and it is con- 
stant for a scale-invariant spectrum. 



In terms of Fourier modes at comoving wavevector k, 
on large scales H(k,t) = TZ pr (k), where ^ pr (k) is the 
constant primordial value of TZ(k), i.e. the value at early 
times sufficiently later than Hubble exit during inflation. 
However, on small scales 1Z decays during radiation domi- 
nation. The total decay incurred through to matter dom- 
ination is described by a linear transfer function T'(fc), 
defined by 



K(k,t m ) =T(fc)ft pr (k), 



(A19) 



where k = |k| and t m is some time during matter domina- 
tion. (T(k) will be distinguished from the temperature T 



APPENDIX B: GAUGE INVARIANCE OF 
ANISOTROPIES 

The calculated temperature anisotropy cannot depend 
on the coordinate choice (in the metric framework) or 
congruence choice (in the covariant framework), since the 
anisotropy is directly observable. A number of workers 
have demonstrated this explicitly in the past in the met- 
ric formalism (see, e.g., [9, 30]). Nevertheless, it will 
be useful to demonstrate the result explicitly using the 
present covariant framework, as it will help to illuminate 
the issues involved. 
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We wish to demonstrate explicitly that the expression 
Eq. (35) for the anisotropy observed at R, namely 



S(H n N)dt + (HStj 



\R ■ 



(Bl) 



is invariant under arbitrary linear gauge transformations, 
if the point R and the four-velocity of the observer are 
held constant. Such a transformation changes the hyper- 
surfaces of constant time according to 



t^t-St, 



(B2) 



for small temporal shift St = St(x fl ). Equivalently, it 
changes the congruence orthogonal to the slicing by the 
spatial gradient of St, 



D»5t, 



(B3) 



at linear order. Under this change in vP, each term in 
the integrand in Eq. (Bl) will change, and the change 
in the slices £j B and implied by Eq. (B2) means 
that the displacements Sto and Sts to the last scattering 
hypersurface Els and to the reception point R will also 
change. 

Straightforward but lengthy calculations using the def- 
initions Eq. (4) to (6) give the following transformations 
under Eq. (B3), to first order: 



0-> 6 + D 2 St + 3HSt, 



(B4) 



1 



a^n^n" -> a^n^n" + n^D^DJt - -D z St, (B5) 



1^ — ► + n^(8i)^, 



(B6) 



where D 2 = D^D^ is the physical Laplacian. Note that, 
by definition, after the gauge transformation (B2) all 
quantities in the integrand in Eq. (Bl) are to be eval- 
uated at the new event temporally displaced by St from 
the original event. This makes no difference at linear or- 
der to quantities that vanish at zcroth order (such as a^ v 

and a M ), but accounts for the term 3HSt in Eq. (B4). 
Next, considering that the quantity Ndt is the proper 
time interval along between hypersurfaces separated 
by coordinate time interval dt, we can easily derive the 
linear transformation law for the lapse perturbation, 



SN -> SN + St. 



(B7) 



Combining Eqs. (B4) to (B7) with Eq. (36), we have the 
linear transformation of the integrand of Eq. (Bl), 



S(H n N) 



S(H n N) + (HSt)+n^n u D tl D l/ 8t 

+ (B8) 



S(H n N) 



dt 



(HSt 



(B9) 



This last line follows from the previous line by straight- 
forward algebra, and contains a coordinate time deriva- 
tive along the null geodesic O. Therefore the integral in 



Eq. (Bl) transforms according to 



S(H n N)dt 



Ie 



S(H n N)dt + (HSt 



(BIO) 

Now all we need are the transformations for the 
boundary terms in the expression for the temperature 
anisotropy, Eq. (Bl). By Eq. (B2), the temporal dis- 
placement (Hdisp between any hypersurface of constant t 
and some fixed hypersurface must transform like 



<^disp — * ^disp — St. 



(Bll) 



Applying this expression to <5tdisp = Str> and <5tdisp = Sts 
gives for the transformation of the boundary terms 



[HSt, 



IR 



(HSti 
(HSt 



<E 
\R ■ 



Combining Eqs. (BIO) and (B12) we finally find 
5T(n») 5T(n^) 



(B12) 



(B13) 



so that the temperature anisotropy is invariant under lin- 
ear transformations of the congruence u^ 1 , or equivalently 
the slicing E t , used in the calculation, if the point of ob- 
servation and the four-velocity of the observer are held 
constant. This of course was to be expected since the 
anisotropics are observable. 

Note that we are free to vary the "background times" 
Ir and Fg at first order. Through Eq. (34) this freedom 
simply shifts the temperature perturbation 5T(n M ) by an 
irrelevant constant. Indeed, we can use this freedom to 
fix 8T(n^) such that its mean over the whole sky vanishes 
for some particular observation point R, 

J ST(n f ")dfl = 0, (B14) 

so that Tr coincides with the mean temperature over the 
sky. Some authors consider it important to make this 
choice (see, e.g., [7]). 



APPENDIX C: GAUGE DEPENDENCE OF 
ANISOTROPIES 

After demonstrating in Appendix B the gauge invari- 
ance of the anisotropies described by Eq. (Bl) for fixed 
observation point R and observer four-velocity, we will 
now show how the anisotropies do depend on the gauge, 
when R and the four-velocity are allowed to transform. 
At linear order, we will see that such transformations 
will only effect the monopole and dipole anisotropies, as 
is well known. 

Consider again the gauge transformation 



t^t-St, 



(CI) 
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with corresponding change in the orthogonal congruence 



D"5t. 



(C2) 



Let us now evaluate the anisotropies using the general 
linear expression, Eq. (Bl), but moving the reception 
point R according to Eq. (CI), and boosting the observer 
four- velocity according to Eq. (C2). Since we move R, 
the corresponding emission point E must also move. If 
R moves to the future, then the corresponding LSS will 
increase in diameter, and E will move radially outwards. 
We can schematically indicate the contributions to the 
change in the anisotropies under Eq. (CI) by 



5E 



d 5T(n») 



dE T R 

d ST(n^) 
+ dvJ I T R 



OR 
5u^. 



6R 



(C3) 



The first term in Eq. (C3) arises due to the change in 
diameter of the LSS (and the entire past light cone). As 
the LSS moves, it samples different perturbation modes, 
so the observed anisotropies change. This effect is great- 
est for the structures at the smallest scales (comoving 
and angular) , and was calculated in detail in [48] . There 
it was shown that this contribution is of order 



d 5T{n^) 
dE T R 



8E 



H{t R )St ST(n^) 
Tr~ 



(C4) 



where 9 is the angular scale of the feature in ques- 
tion. Since St and 5T(n M ) are both first order quantities, 
this effect can be considered to be second order, and so 
will not be considered further here. However, for large 
changes in observation time, substantial changes to the 
anisotropies will be observed [48]. 

To calculate the second and third terms in Eq. (C3), 
note first that by displacing the observation point by St 
and the observation four-velocity by D^St, the displace- 
ments Sto and Sts at R do not change under Eq. (CI): 



6t D (R) -> 5t D (R), 
St B (R) -> St B (R). 



(C5) 
(C6) 



The displacements at the emission point E still transform 
according to Eq. (B12), 



(H8t D + n„St'^)\ E -> (HSt D + n^Si'g) \ 



- (HSt 



(C7) 



Combining Eqs. (C5) to (C7) with the transformation 
Eq. (BIO) for the integral, we find that the anisotropics 
described by Eq. (Bl) transform according to 



ST{n^) 8T{n^) 



Tr, 



(HSt + n^St' 1 *) | 



(C8) 



To illuminate the nature of this change in the 
anisotropics, we can use the multipole expansion of the 
anisotropy, Eq. (38). If we align the polar axis along 
D^St, we have 



HSt 



^HStYooin^), 



n ^ = \/y |^|>io«), 



(C9) 
(CIO) 



where Su^ = D^St. Combining these expressions with 
Eq. (C8), the multipole expansion (38) gives 



a-oo —> a o — V AnHSt, 
/4tt, 

a w — > aio 




(Cll) 
(C12) 



and all other multipolcs arc invariant under the transfor- 
mation. That is, only the monopole and dipole change. 
Therefore the calculation of the higher multipoles is for- 
giving with respect to the care taken regarding gauge. 
However, in Sec. IV, where we calculate the dipole and 
monopole, we must be completely explicit about the spec- 
ification of the frame in which we evaluate the dipole and 
the hypersurface on which we evaluate the monopole. 

To close this discussion of gauge dependence, we note 
that going beyond first order, a boost at the observa- 
tion point R transfers power to all multipoles, and also 
distorts anisotropies through aberration [49] . 
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